Quantum Hall Effect of Dirac Fermions in Graphene: Disorder Effect and Phase 
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We numerically study the interplay of band structure, topological invariant and disorder effect in 
two-dimensional electron system of graphene in a magnetic field. Two distinct quantum Hall effect 
(QHE) regimes exist in the energy band with the unconventional "half-integer" QHE appearing near 
the band center, consistent with the experimental observation. The latter is more robust against 
disorder scattering than the conventional QHE states near the band edges. The phase diagram for 
the unconventional QHE is obtained where the destruction of the Hall plateaus at strong disorder 
is through the float-up of extended levels toward band center and higher plateaus always disappear 
first. We further predict a new insulating phase between f = ±2 QHE states at the band center, 
which may explain the experimentally observed resistance discontinuity near zero gate voltage. 

PACS numbers: 73.43.Cd; 73.40.Hm; 72.10.-d; 72.15.Rn 



With the advances in micromechanical extraction and 
fabrication techniques, high mobility single atomic layer 
of graphite, called graphene, became available recently 1, 
0, M 3j which has attracted much experimental and 
theoretical interest. This new material has many ex- 
traordinary properties such as submicron-scale ballistic 
transport at room temperature, ability to sustain high 
electric currents, and ease of tuning carrier density for 
both electrons and holes. It is generally believed that 
graphene may have potential applications in electronic 
devices 0, |M 0, 01- Graphene is of great fundamen- 
tal interest as well because of its special band struc- 
ture 0. Undoped graphene has one tt electron per lat- 
tice site forming a two-dimensional (2D) extended elec- 
tronic structure. The only states at the Fermi energy 
(Ep — 0) are at two corners of the Brillouin zone, where 
the conduction and valence bands touch. The low en- 
ergy excitations have a linear dispersion relation similar 
to that of the massless Dirac equation, where the elec- 
tron Fermi velocity plays the role of an effective "speed 
of light" . So the electron system of graphene can be 
viewed as a condensed-matter realization of relativistic 
Dirac fermions, and the band touching points are often 
referred to as Dirac points. 

Remarkably, novel quantum Hall effect (QHE) with 
Hall plateaus obeying the unconventional quantization 
rule (Jxy = (fc + l/2)5s \ has recently been observed ex- 
perimentally 0, Q in graphene films in a strong mag- 
netic field. Here, fc is an integer and = 4 stands for 
the spin and sublattice-related degeneracy. In units of 
gs ^ , the "half-integer" quantization of a-j^y was conjec- 
tured to be due to the nontrivial Berry phase of the Dirac 
fermions in a magnetic field jiH,!!. Interestingly cal- 
culations 0, 0, 01 by using an analogy to the 2-1-1 
dimensional Quantum Electro Dynamics have predicted 
a "half-integer" quantized QHE for graphene. Disorder 
effect has been studied by using self-consistent Born ap- 



proximation within the continuous model 12] . However, 
so far the interplay of the band structure, particle-hole 
symmetry, disorder effect, and topological property of 
the energy band in the unconventional QHE has not yet 
been investigated. It is thus highly desirable to perform 
exact numerical calculations by taking into account the 
full band structure and random disorder, in order to re- 
veal the fundamental nature of the new QHE phases and 
related quantum phase transitions in graphene. 



In this Letter, we study the QHE in graphene using a 
tight-binding model. The experimentally observed un- 
conventional Hall plateaus are reproduced around the 
band center Ep = 0, while the conventional integer 
QHE plateaus appear near the band edges. The un- 
usual distribution of the topological invariant quantity- 
Chern number in the energy band and the conservation 
of total Chern number (which is also the geometric Berry 
phase ^3|) for the particle and hole bands account for the 
unconventional "half-integer" QHE. The latter is found 
to be much more stable than the conventional ones near 
the band edges against disorder scattering. We map out 
the whole phase diagram for the QHE and demonstrate 
that the QHE plateaus are destroyed at strong disorder 
(or weak magnetic field) through the float- up of extended 
levels toward the band center. We further identify a new 
insulating phase between v = ±2 QHE states, which 
may explain the experimentally observed discontinuity 
in resistance^, 0| near zero gate voltage. 



We consider a rectangular sample of 2D graphene sheet 
consisting of carbon atoms on a honeycomb lattice [slll^. 
which has totally Ly zigzag chains with atomic sites 
on each chain [lj| . The size of the sample will be denoted 
as Lx X Ly. In the presence of an applied magnetic field 
perpendicular to the graphene plane, the lattice model 
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FIG. 1: (a) Calculated Hall conductance and electron density 
of states in the full energy band for magnetic flux strength </!> = 
||- or M = 48, and (b) the unconventional Hall conductance 
for three different strengths of magnetic flux near the band 
center. Here the disorder strength is set to VF = and the 
sample size is taken to be 96 x 48. 



can be written in the tight-binding form |^ Hsf : 

H = ^ e'"-4cJa+^^«^4c- 



(1) 



where {cia) creates (annihilates) a tt electron of spin 
a on lattice site i with t as the nearest-neighbor hopping 
integral, and Wi is random disorder potential uniformly 
distributed in the interval Wi <E [~W/2, W/2]. The mag- 
netic flux per hexagon cj) = = ^ is proportional 
to the strength of the applied magnetic field B. This 
tight-binding model is valid for describing the full energy 
band and realizes the energy dispersion relation of the 
Dirac fermions near the band center j^. 

The Hall conductance axy can be calculated by using 
the Kubo formula through exact diagonalization of the 
system Hamiltonian fn Fig. la, the Hall conduc- 

tance axy and electron density of states are plotted as 
functions of electron Fermi energy Ep for a clean sam- 
ple W = with magnetic flux = |f , which illustrates 
the overall picture of the QHE in the full energy band. 
According to the behavior of axy, the energy band is nat- 
urally divided into three different regimes. Around the 

2 

band center, axy = v\ is indeed quantized according to 
the unconventional quantization rule v = (k+l/2)gs with 
a degeneracy factor gs = 4 for each Landau level (LL) due 
to two spin directions and two Dirac points. These Hall 
plateaus explain the experimentally observed unconven- 
tional QHE 0, 0] and agree with the results from the 
theory based upon the continuous model |^ 0, 0, ^3 ■ 
In Fig. lb, the Hall conductance in this unconventional 



shown. With decreasing magnetic flux from v — 24 
(f) — more quantized Hall plateaus emerge following 
the same quantization rule. We see that the widths of 
the plateaus are roughly proportional to_v/^ oc V^, in 
agreement with the continuous theory The uncon- 

ventional QHE in the present band model can be under- 
stood in terms of the topological invariant Chern num- 
ber [ilEi. Inside each LL, there are extended states 

characterized by a nonzero Chern integer. The total Hall 
2 

conductance in units of ^ is exactly the sum of the Chern 
numbers of all the occupied extended states mili. The 
additional degeneracy = 4 around the band center 
gives rise to a total Chern number C = 4 for each LL. 

Thus when each additional LL is occupied, the Hall con- 

2 

ductance increments by gs^- At the particle- hole sym- 
metric point Ep = 0, corresponding to the half-filling 
of the central LL, axy — and the total Chern num- 
ber of all the occupied states (hole band) must sum up 
to zero. Now one can count axy from this point, and 
find that the central LL effectively contributes ±(^)^ 
to axy, when Ep is shifted away from the central LL by 

adding particles or holes. This leads to the experimental 

2 

"half-integer" quantization of axy in units of gs^- As 
will be shown below, the total zero Chern number for 
the particle or hole band resulting from the particle-hole 
symmetry in the pure system remains to be true in the 
presence of disorder. 

Near the band edges, each LL carries a total Chern 

number C = 2, and thus the Hall conductance is quan- 

2 

tized as axy — kgs^ with k an integer and gs = 2 for 
spin degeneracy only, which is as same as that in the con- 
ventional QHE systems. Remarkably, around Ep = ±t, 
there are two new critical regions, which separate the 
unconventional and conventional QHE states. The ex- 
tended states in each critical region carry a large nega- 
tive total Chern number, e.g., = for = ||, 
as shown in Fig. la. The existence of the negative Chern 
number regimes around Ep = ±t is crucial for under- 
standing the behavior of the Hall conductance in the 
whole energy band. When the Fermi energy Ep is in- 
creased from the band bottom toward band center con- 
tinuously, following a whole sequence of the conventional 
Hall plateaus, the negative Chern numbers cause a dra- 
matic reduction and a sign inversion of axy, so that the 
unconventional low Hall plateaus with 1^ — —6, —2, 2 ■ ■ ■ 
can reoccur near the band center. This is in contrast to 
the QHE on a square lattice, where stable Hall plateaus 
can only be observed near the band edges 15]. 

We have shown that the unconventionally quantized 
QHE observed in the experiments can be reproduced in 
the lattice model, and is due to the unusual topological 
property of the energy band. We next turn to the effect 
of random disorder on the QHE. In Fig. 2, the Hall con- 
ductance around the band center is shown as a function 
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FIG. 2: Unconventional Hall conductance as a function of 
electron Fermi energy near the band center for four different 
disorder strengths each averaged over 200 disorder configura- 
tions. Inset: conventional Hall conductance near the lower 
band edge. Here M — 96 and the sample size is 96 x 48. 

of Ep for magnetic flux = |f and four different values 
of disorder strength W. We see that with increasing W, 
higher Hall plateaus (with larger are destroyed first. 
At W — O.bt, the plateaus with v — ±10, ±6 and ±2 re- 
main well quantized, while at = 2.Qt all the plateaus 
except for the v — ±2 ones are destroyed. These last 
two plateaus will eventually disappear at W ^ 2.5t. For 
comparison, the QHE near the lower band edge is shown 
in the inset, where all plateaus are found to be destroyed 
at a much weaker disorder W = l.Qt. This clearly in- 
dicates that, under the same conditions, the unconven- 
tional QHE is much more stable than the conventional 
one. This is attributed to the Dirac-fermion-like linear 
dispersion relation around the band center, where the 
widths of the LL gaps are proportional to \/B instead of 
B. We also notice that, axy always vanishes at Ep = 
for all W, due to the fact that the whole particle or hole 
band carries zero total Chern number as a topological 
invariant in the disordered system. 

We further study the quantum phase transition of the 
graphene electron system and establish the phase dia- 
gram for the QHE. This can be done relatively conve- 
niently by calculation of the finite-size localization length 
A on an essentially infinitely long bar of width Ly (length 
Lx > 10^) by using the well-established recursive Green's 
function approach • We present the calculated phase 
diagram in Fig. 3a, for a relatively large magnetic flux 
(/) = ll for clarity, while the topology of the phase di- 
agram is essentially universal, independent of 0. In the 
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W — Ep plane, different QHE plateaus with axy — v\ 
are separated by extended states, where A grows linearly 
with increasing bar width Ly. With the increase of W , 
each plateau can be destroyed through a transition v ^ Q 
to the insulating phase and higher plateaus disappear 
first. In Fig. 3b-3d, we show examples of calculated lo- 
calization length to explain how the phase boundaries 
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FIG. 3: (a) Phase diagram for the unconventional QHE 
regime in graphene at M = 48, which is symmetric about 
Ef ~ 0. (b) to (d): Normalized localization lengths calcu- 
lated for three bar widths Ly — 48, 96 and 144, as the phase 
boundary is crossed by the paths indicated by the arrows B, 
C and D in (a), respectively. 



in the phase diagram are determined. In Fig. 3b, the 
normalized localization length \/Ly for Ep = —0.3/: and 
three sample widths Ly = 48, 96 and 144 is plotted as a 
function of W, which corresponds to a v = —2 tran- 
sition as indicated by the arrow B in the phase diagram. 
The sample length Lx ranges from 10® up to 5 x 10®, 
so that the relative error due to statistical fiuctuations 
in A reduces to about 2%. We see clearly that X/Ly 
is peaked a,t W = Wc — 2.7t, an indication of an ex- 
tended critical point separating the ly = —2 plateau from 
the outside insulating phase. Finite-size scaling con- 
firms that the localization length at the thermodynamic 
limit becomes divergent at Wc- This is consistent with a 
new "float-up" picture [isj . where some negative Chern 
number states are coming from lower energy, and moving 
toward the band center with increasing W, which sweep 
across Ep = —0.3/ at W = Wc, causing the collapse of 
the v = —2 plateau. 

Figure 3c shows the normalized localization length as 
a function of W at Ep — —0.55/, corresponding to the 
path indicated by the arrow C in the phase diagram Fig. 
3a. We see that a peak of localization length occurs at 
W = Wc — 1.6/. Similarly to Fig. 3b, the peak indi- 
cates the destruction of the — —6 QHE state and its 
transition into the insulating phase. However, we note 
that here the localization length A is relatively large, be- 
ing much greater than the largest Ly that is reachable in 
our calculations. So X/Ly does not decrease visibly with 
increasing Ly, and one cannot rule out a possibility that 
the higher plateau to insulator transition happens in a 
critical region with a small finite width AWc instead of 



4 




0.5 I ' I ' ' ' I ' ' ' ^ ' ' ' ^ ' ' ' ^ ' L 

0.8 1.6 2.4 3.2 4.0 
W/t 



FIG. 4: The normalized localization length X/Ly for three 
bar widths Ly — 48, 96 and 144, when the phase boundary is 
crossed with varying disorder strength W at Ef ~ — O.Olt, as 
indicated by the arrow E in Fig. 3a. Inset shows the localiza- 
tion length ^ at the thermodynamic limit determined through 
one-parameter scaling. 

at a critical point Wc, where electron states become de- 
localized. All the phase boundaries separating the QHE 
phases from the insulating phase at strong W, indicated 
by the solid line with open circles in Fig. 3a, are deter- 
mined in the same way. 

To determine the phase boundary between different 
QHE states, the localization length is calculated as a 
function of electron Fermi energy Ep for fixed W. As 
shown in Fig. 3d, corresponding to the path indicated 
by the arrow D in Fig. 3a, a peak in X/Ly occurs at 
Ep = Ec — — 0.46i, which indicates a critical point sepa- 
rating v — —6 and —2 plateaus. All the phase boundaries 
indicated by dotted lines with cross symbols are deter- 
mined in the same manner. 

An important feature of the phase diagram is that the 
v = —2 and 2 plateaus around the band center are no 
longer connected to each other in the presence of dis- 
order, separated by a new insulating phase in between. 
Corresponding to such a phase transition, along the path 
indicated by the arrow E in Fig. 3a, we show X/Ly at 
Ep ~ — O.Oli as a function of W in Fig. 4. Clearly 
a largely reduced critical disorder strength Wc — l.Ot 
is observed in Fig. 4. Following the standard finite-size 
scaling analysis, we find that all the data in Fig. 4 can 
be well fitted by a one-parameter scaling relation ^3 
X/Ly = f{Ly/Cl for Ly = 48,96 and 144. The fitting 
parameter £,{W) is the localization length at the thermo- 
dynamic limit. Its value is plotted in the inset of Fig. 4 as 
a function of W , which becomes divergent at Wc — \M. 
Interestingly, the splitting indicates a singularity in the 
resistance near Ep — 0, as an insulating phase is char- 
acterized by a divergent resistance pxx, while a plateau 
state has zero or finite pxx (being finite only at the criti- 
cal point) at low temperature limit. This feature may ex- 
plain the divergent trend in pxx and discontinuity in pxy 



observed near zero gate voltage in the experiments [a, 13 • 
We have also confirmed the above phase boundaries 
by calculating the Thouless number 18], which is pro- 
portional to the longitudinal conductance. In particular, 
we observed that the Thouless number as a function of 
Ep at fixed W shows two peaks near the band center, 
while a dip occurs at Ep = 0, which is consistent with 
the splitting of the = — 2 to 2 transition with a new 
insulating phase emerging near Ep = 0. 

In summary, we have numerically investigated the 
QHE in 2D graphene based upon a lattice model. The ex- 
perimentally discovered unconventional quantization of 
QHE is reproduced near the band center, which is under- 
stood in terms of the novel distribution of the topological 
Chern integers in the energy band. The phase diagram 
indicates a new float-up picture, in which the extended 
levels move toward band center with increasing disorder 
strength, causing higher plateaus to disappear first. The 
unconventional QHE plateaus around the band center 
are found to be much more stable than the conventional 
ones near the band edges. A new insulating phase is pre- 
dicted to emerge at the band center, between two v = ±2 
QHE states, which is consistent with the experimentally 
observed resistance discontinuity near zero gate voltage. 

Acknowledgment: This work is supported by ACS- 
PRF 41752-AClO, Research Corporation Fund CC5643, 
the NSF grant/DMR-0307170 (DNS), a grant from the 
Robert A. Welch Foundation under the grant no. E- 
1146 (LS), and the NSFC grants 10374058 and 90403016 
(ZYW). 



[1] K. S. Novoselov et al., Science 306, 666 (2004). 

[2] C. Berger et al., 3. Phys. Chem. B 108, 19912 (2004). 

[3] Y. Zhang, J. P. SmaU, W. V. Pontius and P. Kim, Appl. 
Phys. Lett. 86, 073104 (2005); Y. Zhang, J. P. SmaU, 
E. S. Amori and P. Kim, Phys. Rev. Lett. 94, 176803 
(2005). 

[4] J. S. Bunch et al., Nano Lett. 5, 287 (2005). 

[5] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

[6] K. S. Novoselov et al.. Nature 438, 197 (2005). 

[7] Y. Zhang, Y.-W. Tan, H. L. Stormer, and Philip Kim, 

Nature 438, 201 (2005). 
[8] M. Wilson, Physics Today 59, 21 (2006). 
[9] V. P. Gusynin and S. G. Sharapov Phys. Rev. Lett. 95, 

146801 (2005). 

[10] N. M. R. Peres, F. Guinea and A. H. C. Neto, 

cond-mat/0506709 
[11] E. McCann and V. I. Fal'ko, cond-mat/0510237 
[12] Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 
[13] D. J. Thouless, et. al, Phys. Rev. Lett. 49, 405 (1982). 
[14] L. Sheng, D. N. Sheng, C. S. Ting and F. D. M. Haldane, 

Phys. Rev. Lett. 95, 136602 (2005). 
[15] D. N. Sheng and Z. Y. Weng, Phys. Rev. Lett. 78, 318 

(1997); D. N. Sheng, Z. Y. Weng and X. G. Wen, Phys. 

Rev. B 64, 165317 (2001). 
[16] Y. Huo and R. N. Bhatt, Phys. Rev. Lett. 68, 1375 



5 



(1992). [18] J. T. Edwards and D. J. Thouless, J. Phys. C 5, 807 

[17] A. MacKinnon and B. Kramer, Phys. Rev. Lett. 47, 1546 (1972). 
(1981); Z. Phys. B 53, 1 (1983). 



